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ABSTRACT 

We report the first extensive study of stellar Rb abundances. High-resolution 
spectra have been used to determine, or set upper limits on, the abundances of this 
heavy element and the associated elements Y, Zr, and Ba in 44 dwarfs and giants 
with metallicities spanning the range —2.0 < [Fe/H] < 0.0. In metal-deficient stars 
Rb is systematically overabundant relative to Fe; we find an average [Rb/Fe] of +0.21 
for the 32 stars with [Fe/H] < — 0.5 and measured Rb. This behavior contrasts with 
that of Y, Zr, and Ba, which, with the exception of three new CH stars (HD 23439A 
and B and BD+5°3640), are consistently slightly deficient relative to Fe in the same 
stars; excluding the three CH stars, we find the stars with [Fe/H] < —0.5 have average 
[Y/Fe], [Zr/Fe], and [Ba/Fe] of -0.19 (24 stars), -0.12 (28 stars), and -0.06 (29 stars), 
respectively. The different behavior of Rb on the one hand and Y, Zr, and Ba on 
the other can be attributed in part to the fact that in the Sun and in these stars Rb 
has a large r-process component while Y, Zr, and Ba are mostly s-process elements 
with only small r-process components. In addition, the Rb s-process abundance is 
dependent on the neutron density at the s-processing site. Published observations of 
Rb in s-process enriched red giants indicate a higher neutron density in the metal-poor 
giants. These observations imply a higher s-process abundance for Rb in metal-poor 
stars. The calculated combination of the Rb r-process abundance, as estimated for 
the stellar Eu abundances, and the s-process abundance as estimated for red giants 
accounts satisfactorily for the observed run of [Rb/Fe] with [Fe/H]. 

Subject headings: stars: abundances — stars: Population II — Galaxy: halo — nuclear 
reactions, nucleosynthesis, abundances 

1. INTRODUCTION 

In the solar neighborhood reside stars of differing metallicity. Stars with a metallicity 
of approximately the solar value belong overwhelmingly to the Galactic disk. Stars of lower 
metallicity are on orbits that identify them as members of the Galactic halo. The metallicity 
[Fe/H] ~ — 1 may be taken as a rough boundary between disk and halo stars. One tool for 
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unravelling the evolution of the Galaxy is the measurement of the chemical compositions of stars 
as a function of metallicity. In this paper, wc present the first extensive series of measurements of 
the rubidium abundance in disk and halo stars. Rubidium is potentially a special diagnostic of the 
neutron capture s-process. 

Rubidium is present in two isotopic forms: ^^Rb which is stable and ^^Rb which with a 
half-life of 5 x 10^'^ yr may be deemed effectively stable from the astrophysical point of view. As 
we remark later, astronomical detection of Kb must rely on the Rb l lines that do not permit 
measurement of the relative isotopic Rb abundances from stellar spectra. Hence, we discuss the 
elemental Rb abundance. Analysis of the solar system abundances of Rb and adjoining elements 
shows that the neutron capture s- and r-processes are about equally responsible for the synthesis of 
Rb. Scrutiny of the s-process abundances shows that the 'main' s-process not the 'weak' s-process 
is the principal source of Rb's s-process component. The main s-process which manufactured 
elements heavier than about Rb is identified with the He-burning shell of intermediate and low 
mass AGE stars. The weak s-process is identified with He-core and possibly C-core burning of 
massive stars. 

Evolution of the Galaxy's s and r-process products is rather directly observed from the stellar 
abundances of elements that are predominantly attributable to either the s or to the r-process. 
Traditional tracers include Ba for the s-process and Eu for the r-process. Elucidation of the 
operation of the neutron capture processes requires observations of more than a single element 
per process. As an example, we note that the abundance ratio of a 'light' to a 'heavy' s-process 
elements, say Zr to Ba, provides information on the integrated exposure of material to neutrons. 
Rubidium with a roughly equal mix of s and r-process contributions, and an unfavorable electronic 
structure for ready detection in stellar spectra would seem to be an element of little interest. 
Closer inspection of the working of the s-process shows, however, that Rb offers a unique insight 
into the process: Rb's role as a monitor of the neutron density at the s-process site. 

Along the s-process path, Rb is preceded by krypton with the path entering Kr at ^''Kr and 
exiting at either ^^Kr or ^*^Kr. Unstable ^^Kr controls the exit. At low neutron density at the 
s-process site, stable '^'^Kr is converted by neutron capture to ^^Kr that decays to stable ^^Rb with 
the path continuing to ^^Sr. At high neutron densities, ^^Kr does not /3-decay but experiences 
a neutron capture and is converted to stable ^^Kr. Subsequent neutron capture by ^^Kr leads 
by /3-decay of ^^Kr to (effectively) stable ^^Rb. When a steady flow along the s-process path is 
attained, the density of a nuclide is given approximately by the condition that djiVj ~constant 
where cTj is the neutron capture cross-section of nuclide i and Ni is the abundance of that nuclide. 
Since as7 — agb/W for the Rb isotopes, the switch of the ^^Kr branch from its low neutron density 
routing through ^^Rb to its high neutron density routing through ^^Rb increases the total Rb 
abundance by about an order of magnitude relative to the abundance of other elements in this 
section of the s-process path, such as Sr, Y, and Zr. The isotopic mix of Rb is obviously altered as 
a function of neutron density but this is not measureable for cool stars. (Krypton is undetectable 
spectroscopically in cool stars.) Operation of the ^^Kr branch is more complicated than sketched, 



-3- 



for example, neutron capture from ^^Kr feeds not only the ^^Kr ground state but a short-lived 
isomeric state that at all reasonable neutron densities provides some leakage to *^Rb. A thorough 
discussion of the ^"''Kr branch is provided by Beer & Macklin (1989) and its use in determining the 
effective neutron density of the s-process in stars is discussed by Tomkin & Lambert (1983) and 
Lambert et al. (1995). When detailed abundance measurements are available as in the case of the 
carbonaceous chondrites, several branches along the s-process path serve as neutron densitometers 
but Rb is the only low neutron density branch available to stellar spectroscopists. (At higher 
neutron densities, a branch controlled by ^^Zr is exploitable in cool stars showing ZrO bands 
[Lambert et al. 1995].) 

A primary reason for the near neglect of Rb in reports on quantitative spectroscopy of stars 
is that it is a trace alkali element. The Rb atom's low ionization potential (4.177 eV) ensures 
that Rb is primarily ionized but the rare-gas electronic structure of Rb+ provides resonance lines 
in the far ultraviolet. Detection of Rb via the Rbi resonance lines at 7800 and 7947 A at the 
expected low Rb abundances is possible for cool dwarfs and giants, as our exploratory synthetic 
spectra indicated. Stars for observation at high-spectral resolution were selected from Schuster & 
Nissen's (1988) catalog of dwarfs, and from Pilachowski, Sneden, & Kraft's (1996) list of giants. 
Emphasis was placed on metal-poor stars such that the metallicity range —2 < [Fe/H] < —0.5 is 
well represented but metallicities [Fe/H] > —0.5 are poorly represented. The following sections 
describe the observations, the method of analysis, the results, and present an interpretation of the 
Rb abundances relative to the abundances of other elements (Fe, Y, Zr, and Ba) obtained here. 

2. OBSERVATIONS AND DATA REDUCTION 

The program stars are listed in Table 1. They comprise 32 G and K dwarfs and subgiants with 
metalhcitics of -1.8 < [Fe/H] < 0.0 and 12 G and K giants with metallicities of -2.0 < [Fe/H] 
< —0.6. The observations were made at the McDonald Observatory with the 2.7-m telescope and 
2dcoude echelle spectrometer (Tull et al. 1995). All the program stars were observed at the F3 
focus at a resolving power of R = 60,000 and, in addition, eight of the brighter stars were also 
observed at the Fl focus at a resolving power of 200,000. In order to minimise the influence of 
cosmic rays, two observations in succession, rather than one longer observation, were generally 
made of each star. Four different detectors were used for the observations with resolving power 
60,000: a Texas Instruments CCD with 15 ;Um^ pixels in a 800x800 format, a Tektronix CCD 
with 27/ini^ pixels in a 512x512 format, the Goddard Advanced Imaging System CCD with 
21 /xm^ pixels in a 2048x2048 format, and a Tektronix CCD with 24 ^um^ pixels in a 2048x2048 
format. The first two of these CCDs provided partial coverage of the wavelength interval ~ 5500 
- ~ 8000 A with large gaps between the end of one spectral order and the beginning of the next. 
The last two CCDs, which are much larger, provided nearly complete coverage from 4000 - 
~ 9000 A; coverage was complete from the start of this interval to 5600 A and substantial, but 
incomplete, from 5600 A to the end of the interval. The Tektronix CCD with 24 /xm^ pixels 
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in a 2048x2048 format was used for the 200,000 resolving-power observations, which provided 
partial coverage of the region from ~ 5500 - ~ 8000 A. The typical signal-to-noise ratio of the 
extracted one-dimensional spectra is between 100 and 300 at red and near-infrared wavelengths 
for the 60,000 resolving-power observations, while it is typically between 100 and 250 at the same 
wavelengths for the 200,000 resolving-power observations. 

The only accessible Rb lines in stellar spectra are the two Rbi resonance lines at 7800.3 
and 7947.6 A. Since these lines are typically weak we concentrated our attention on the stronger 
7800.3 A line and did not pursue the 7947.6 A line, which is half as strong as the 7800.3 A line. 
Figure |l] shows examples of the 7800.3 A Rb i line in the program stars; as may be seen, the line is 
partially blended with a stronger Sii line at 7800.0 A. 

The data were processed and wavelength calibrated in a conventional manner with the IRAF]] 
package of programs on a SPARC5 workstation. The 7800.3 A Rbi line was analyzed by spectrum 
synthesis because of the presence of the Sii line. Equivalent widths were measured for lines of 
other elements used in the investigation; these were lines of Fe i and Fe ii and the available lines 
(Yii, Zri, Ball, and Ndii) of other heavy elements besides Rb. Lines suitable for measurement 
were chosen for clean profiles, as judged by inspection of the solar spectrum at high resolution and 
signal-to- noise ratio (Kurucz et al. 1984), that could be reliably measured in all, or most of, the 
program stars. Moore, Minnaert, & Houtgast (1966) was our primary source of line identification. 
The equivalent width of each line was measured with the IRAF measurement option most suited to 
the situation of the line; usually this was the fitting of a single, or multiple, Gaussian profile to the 
line profile. These equivalent widths were measured from the spectra with 60,000 resolving power; 
the 200,000 resolving-power spectra were not used because their much more limited wavelength 
coverage excluded most lines of interest. Table 2 gives basic information for the Rbi line and the 
lines of the other elements. The equivalent widths of the lines are available at JT's World Wide 
Web site ( http :/ /anchor . as . ut exas . edu/ 1 omkin . html] ) . 



The spectrum of an asteriod (Iris), observed with the same instrumental setup as that used 
for the 60,000 resolving-power observations of the program stars and reduced and measured in the 
same manner, provided solar equivalent widths for these lines. 



3. ANALYSIS 

An LTE model-atmosphere abundance analysis was made relative to the Sun. Here we briefly 
discuss the selection of atomic data for the lines, the model atmospheres, and the abundance 
determinations for Rb and the other heavy elements. 



^ IRAF is distributed by the National Optical Astronomical Observatories, which is operated by the Association 
for Universities for Research in Astronomy, Inc., under contract to the National Science Foundation. 
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3.1. Line Data 

The lines used for the abundance determinations are given in Table 2. Our first choice for 
^/-values for the lines was modern laboratory ^/-values; basic data for the lines, including the 
5/- values and their sources, are given in Table 2. In particular we note that the 5/- value of the 
7800.3 A Rbl Hue, log 5/ = +0.13±0.04 (Wiesc k Martin 1980), is rchably determined. For some 
lines, for which reliable ^/-values are not available, wc used solar (^/-values instead. These were 
calculated with the solar equivalent widths given in Table 2, the solar atmosphere of Holwegcr & 
Miiller (1974), a microturbulence of 1.15 km s~^ (Tomkin et al. 1997), and the abundances given 
in Table 2. The solar equivalent widths were measured from the spectrum of the asteriod Iris. 

The adoption of the laboratory ^/-values for the Fe I lines follows the prescription of Lambert 
et al. (1996). In particular, we make a small correction (see Table 2) to the 5/- values of May, 
Richter, &; Wichelmann (1974) to normalize them to those of O'Brian et al. (1991) and Bard, 
Kock, &; Kock (1991). No correction is needed to put the solar Fel line ^/-values on the same 
scale as the laboratory ^/-values; we find: log (//(solar) - log (//(lab) = +0.03ib0.05 (4 lines, 6 
(//-values). Accurate experimental (//-values are not available for our Fell lines, therefore we have 
used solar (//-values for these lines. 

Although the three Zrl lines all have modern laboratory (//-values (Biemont et al. 1981), 
we find their solar ^'/-values are significantly larger than the laboratory ones: log y/(solar) - log 
5'/(lab) = -|-0.41±0.04. Our adopted solar Zr abundance — log e(Zr) = 2.60 (Anders & Grevesse 
1989) — is not a factor in this discrepancy because it is based on Biemont et al.'s 5/- values and is 
very similar to the value — log e(Zr) = 2.56 — Biemont et al. derived in their own investigation of 
the solar Zr abundance. Some of the discrepancy is attributable to the use of different ionization 
potentials for Zrl; we have used the accurate value of 6.634 eV (Hackett et al. 1986) while Biemont 
et al. must have used the significantly higher old value of 6.84 (Allen 1973). This accounts for 
0.20 dex of the discrepancy. Line-to-line variation of individual line abundances may also account 
for some of the discrepancy; the individual line abundances Biemont et al. derive for these three 
lines are on average 0.06 dex larger than the average Zr abundance they determine from the Zr I 
lines. Although the remaining 0.15 dex discrepancy is not readily accounted for, the manner in 
which these three lines strengthen in the cooler program stars leaves no doubt that they are low 
excitation lines of a neutral species and is thus consistent with their identification as Zr i lines. In 
order to minimize the influence of whatever is causing the residual discrepancy between the solar 
and laboratory ^/-values for these lines, we have chosen to adopt the solar ^/-values for them in 
our analysis. 

Two Yll lines (5200.42 and 5402.78 A) were rejected when it became evident, during an 
initial abundance analysis of all the stars, that they gave significantly larger Y abundances than 
the other Yll lines in many of the stars. 

Atomic data and solar equivalent widths for the lines are given in Table 2. 



-6- 



3.2. Model Atmospheres 

Plane parallel, line blanketed, flux constant, LTE MARCS model atmospheres, which derive 
from a development of the programs of Gustafsson et al. (1975), were used for the abundance 
analysis. 

The determination of parameters (effective temperature, surface gravity, metallicity, and 
microturbulence) for the model atmospheres was done in two steps. First we chose preliminary 
parameters for each star which were used to calculate an initial set of model atmospheres. We 
then used the initial model atmospheres and the equivalent widths of the Fe I and Fe II lines to 

iteratively adjust the atmospheric parameters to determine an adopted set of parameters and 
atmospheres that are consistent with the Fei and Fell line data. We now briefly describe these 
two steps. 



3.2.1. Choice of Preliminary Atmospheric Parameters 

Stromgren photometry provided the primary means of determining preliminary parameters 
for the dwarfs and subgiants, which are mostly from the catalogue of Schuster & Nissen (1988) 
supplemented by three stars from Carney et al. (1994) and four stars from the Bright Star 

Catalogue (Hoffleit & Jaschek 1982). The stars are intrinsically faint and nearby; an examination 
of their parallaxes (see below) shows they are all within 100 pc, except for two stars which are at 
llliblT and 137±18pc. The interstellar reddening of the stars is therefore negligible (Schuster &; 
Nissen 1989a) so we have not corrected their indices. 

Our chief source of preliminary effective temperatures for the dwarfs and subgiants was the 
color index b — y coupled with the Tgff vs. b — y calibration of Alonso, Arribas, & Martmez-Roger 
(1996a, equation 9), who used the infrared flux method to determine effective temperatures for 
metal-deficient stars. For the components of the visual binary HD 23439, which do not have 
their own b — y, and for HD 150281, which also does not have b — y, we took the effective 
temperatures from Carney et al. (1994); for 61 Cyg A and B, which are too cool for the applicable 
range of Alonso et al.'s calibration, we took the effective temperatures from Alonso, Arribas, &: 
Martmez-Roger (1996b). 

Initial surface gravities for the dwarfs and subgiants were determined by the relations 
g oc M/R'^ and L oc R'^T^q^ where M is the stellar mass, R the radius, and L the luminosity, with 
the luminosities being set by the Hipparcos parallaxes. We followed the prescriptions of Nissen, 
H0g, & Schuster (1997), who have successfully applied this method to determine surface gravities 
for 54 metal-poor stars. Most of the stars have Hipparcos parallaxes and all of these parallaxes are 
of sufficient accuracy; the largest uncertainty (c/vr) is 0.15 with most uncertainties being much 
smaller than this. For the small number (five) of stars without Hipparcos parallaxes we took 
trigonometric parallaxes from Gliese &; Jahreiss (1979) or determined photometric parallaxes from 
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V and My, with My estimated from the uvby photometry and the recipes of Nissen & Schuster 
(1991). 

PreUminary metaUicities for the dwarfs and subgiants were estimated from the uvby — /? 

photometry and the caHbration of Schuster & Nissen (1989b, equation 3). For a small number 
(five) of stars the results of Carney et al. (1994) or Alonso et al. (1996b) were used instead. 

An initial microturbulence of 1.0 km s""^, which is representative of dwarfs in this temperature 
range (Feltzing &; Gustafsson 1998), was used for the dwarfs and subgiants. 

For the giants, which arc taken from Pilachowski et al.'s (1996) medium-resolution 
spectroscopic investigation of Na abundances in metal-poor giants, we adopted Pilachowski et 
al.'s atmospheric parameters as initial parameters. Pilachowski et al.'s effective temperatures 
and gravities are photometrically based as modified by their spectroscopic results, while their 
metaUicities and microturbulences are from their spectroscopic analysis. 

3.2.2. Determination of Adopted Atmospheric Parameters 

Model atmospheres were computed using the MARCS code (Gustafsson et al. 1975). Those 
for the dwarfs and subgiants were calculated by interpolation in a grid of MARCS models, which 

spanned the range of dwarf and subgiant parameters and was provided by B. Edvardsson. The 
models of Pilachowski et al. (1996), who also used MARCS models for their abundance analysis, 
were used as preliminary models for the giants and were provided by C. Sneden; iterations of the 
giant models with modified atmospheric parameters were calculated directly using the MARCS 
code. 

Models with the preliminary parameters and the line analysis code MOOG (Sneden 1973) were 
then applied to the equivalent width data for the Fei and Fell lines. Trends of the Fei abundances 
with line excitation potential were used to check the preliminary effective temperatures and trends 
of the Fe l abundances with equivalent width were used to check the preliminary microturbulences. 
Where necessary revised parameters were determined, new models calculated, and a new round of 
abundance calculations done for the Fei and Fell lines. Next the Fc abundances from this round 
of calculations were used to calculate a new set of models and do another round of abundance 
calculations for the Fe lines. In a final round of calculations the surface gravities were adjusted, 
and new models were calculated, so the Fei and Fell lines gave the same Fe abundance. These 
final adopted parameters for the program stars are given in Table 1. 

The adopted parameters are generally only moderately different from the initial parameters. 
For the dwarfs and subgiants the analysis of the Fe lines led to revised effective temperatures for 
13 stars and an average temperature increase for these stars of 140ib75 K (d of the individual 
differences). The revised gravities of the dwarfs and subgiants tend to be lower than the 
preliminary gravities, but the inconsistency is small; for the 28 dwarfs and subgiants with 
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measurements of both Fe i and Fe ii lines, for which spectroscopic gravities can thus be determined, 
the average downward revision of the log g is 0.12ib0.12 [a of the individual differences). This 
suggests that there is no serious inconsistency between the gravities of these stars and the four 
remaining dwarfs and subgiants for which we adopt preliminary gravities because they have no 
Fell lines. 

It is of interest to see how the differences that we find between our spectroscopic gravities and 
the preliminary Hipparcos-based gravities compare with what Allende Prieto et al. (1999) found 
in a thorough examination, of spectroscopic gravities, for nearby stars, taken from the literature 
versus Hipparcos-based gravities. We confine the comparison of log g^pcc ^ log ^Hipp for our results 
and theirs to the temperature range (4900 - 5500 K) of our dwarfs and subgiants. (This range 
excludes the four coolest dwarfs and subgiants because they do not have any measured Fell lines.) 
As mentioned earlier, we find an average difference log ^spec ~ log 5Hipp = -0.12±0.12 (28 stars, a 
of individual differences) for our dwarfs and subgiants, while Allende Prieto et al. find an average 
difference of ~0.26=b0.29 (9 stars) for their sample of stars. Our results and Allende Prieto et al.'s 
thus both show that the spectroscopic gravities tend to be smaller than the Hipparcos gravities. 
Also, the consistency between the spectroscopic and Hipparcos gravities of our stars is somewhat 
better than it is for Allende Prieto et al.'s sample of stars.0 

Revisions of effective temperature were required for five of the giants; the average temperature 
increase for these stars was 100ibl40K. All the giants have measurements of both Fei and Fell 
lines thus allowing revision of their preliminary gravities for all 12 stars; the average change of 
log g was -0.04zb0.37. The adopted effective temperatures and gravities of the giants thus are in 
good agreement with the preliminary values (Pilachowski et al. 1996). Our results, which use 
higher resolution spectra and more numerous Fe lines than those of Pilachowski et al., thus confirm 
their results. Our Fe abundances are also in good agreement; the average difference between our 
[Fe/H] determinations and theirs is -0.07=b0.10. 

We conclude this section with a brief discussion of the potential influence of non-LTE on 
the effective temperatures derived from the Fe i line excitation. As remarked earlier, the primary 
source of our initial effective temperatures (for the dwarfs and subgiants) was b — y and the 
Teff vs. b — y calibration of Alonso, Arribas, & Martinez- Roger (1996a, ), where the effective 
temperatures in their calibration were determined by the infrared flux method. The b — y based 
effective temperatures are thus free of non-LTE effects. One way to estimate the possible influence 
of non-LTE on the effective temperatures derived from the excitation of Fei lines, therefore, 
is to consider the difference between the excitation-based temperatures and the b — y-based 
temperatures. Of course non-LTE effects are not the only possible source of such a difference so 
this check is indicative, rather than conclusive. The average difference T^g (Fe i) — Tefj {b — y) 



^ We have not determined Hipparcos-based gravities for the giants because, with the exception of Arcturus, they 
are much more remote than the dwarfs and subgiants. Although they aU have Hipparcos parallaxes, the errors in the 
parallaxes are comparable to the parallaxes for most of them. 
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= +45±68K (cr of the individual differences, 20 stars), where the calculation includes not only 
stars whose initial effective temperatures were revised, but also stars for which no revision was 
necessary - as long as they had enough Fcl lines to define the line excitation. This difference is 
small and indicates that any non-LTE influence on the determination of temperatures from the 
Fe I lines is minor. 

3.3. Abundance Determinations for Rb and the Other Heavy Elements 

Abundances were determined by matching the observed line strengths and theoretical line 
strengths calculated by MOOG with the adopted model atmospheres. As remarked earlier, the 
Rbl line was treated by means of spectrum synthesis, while the lines (Yll, Zrl, Ball, and Ndll) 
of the other heavy elements were treated by means of equivalent widths. 

The spectrum synthesis of the 7800.29 A Rbl line includes the hyperfine structure of the 

®^Rb and ^^Rb isotopes, each of which is split into two components, and the blending Si I line at 
7800.00 A. The accurately known wavelengths and relative line strengths of the hyperfine structure 
components were taken from Lambert &: Luck's (1976) analysis of Rb in Arcturus. We adopted a 
terrestrial abundance ratio (^^Rb/^^Rb = 3) for the Rb isotopes. Although in principle it would 
be desirable to make direct measurements of the stellar isotopic Rb abundances from the exact 
shape of the Rbl line profile, in practice extreme departures from the terrestrial isotope ratio are 
required before there is appreciable distortion of the line profile. Lambert & Luck, for example, 
found that for their spectra, which had a resolving power of 195,000 that is similar to the high 
resolution spectra of the present investigation, the ®^Rb/®^Rb ratio had to be as low as 1 or as 
high as 10 to cause even a small variation of the line profile; they concluded that the isotope 
ratio in Arcturus is terrestrial with a large uncertainty. Direct determination of Rb isotopic 
abundances, therefore, is not the thrust of our investigation, although we note that comparison 
of the observed and synthesised spectra of the Rb line in our program stars does not show any 
variations attributable to non-terrestrial isotope ratios. We also note that the indeterminacy of 
the Rb isotope ratios does not interfere with the measurement of the elemental Rb abundances; 
the Rb I line is weak in all the program stars so the Rb abundances it provides are not affected by 
changes of the isotopic mixture. 

The synthesis of the Rbl line also includes the Si I line to the blue. No reliable experimental 
oscillator strength is available for the Si I line so a solar oscillator strength (log gf = -0.65) was 
adopted. The instrumental and macroturbulent broadening, as well as thermal and microturbulent 
broadening, were included in the synthesis. The macroturbulent broadening was set by matching 
the profile of the clean nearby Ni I line at 7797.6 A. 

Spectrum synthesis of the solar Rbl line, using the solar model of Holweger & Miiller (1974) 
and the Kurucz et al. (1984) solar atlas, provides a solar Rb abundance log e(Rb) = 2.60ib0.07, 
which is the same as the Rb abundance in Anders & Grevesse's (1989) compilation of solar 
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abundances. We note, however, the discrepancy between the photosperic abundance and the 
somewhat lower meteoritic abundance of 2.40ib0.03 (Anders & Grevesse 1989). Although it might 
be speculated that the relatively low melting and boiling points of Rb (39 and 688 C, respectively) 
may make it behave like a volatile element and so explain the low meteoritic abundance, this does 
not appear to be the case. Potassium, which is isoelectronic with Rb and has only slightly higher 
melting and boiling points (63 and 759 C, respectively), shows no discrepancy of its photospheric 
and meteoritic abundances, which are log e(K) = 5.12ib0.13 and 5.13ib0.03, respectively. Although 
the discrepancy is a potential source of concern, we note that there is no evidence of the Rb i line 
being affected by an unknown blend; in particular, in our program stars the strengthening of the 
line with decreasing effective temperature is consistent with the behavior of a resonance line of a 
heavy-element neutral species. 

The observed and synthesised Rb i line profiles for a sample of stars are shown in Figure ||. 
Table 3 gives the abundances of Rb and the other heavy elements. 

The Rb abundances derived from observations made at the Fl focus (resolving power = 
200,000) and the F3 focus (resolving power = 60,000) are highly consistent; for stars observed at 
both foci the average difference between the Fl- and F3-based [Rb/H] is -0.03ib0.01 (s.e., five 
stars). For bright stars, such as Arcturus and fi Gas, the greater detail provided by the higher 
resolution Fl observations allows for more precise determination of the Rb abundance. In fainter 
stars there is not much to choose between the Fl and F3 observations because the greater spectral 
detail of the Fl observations tends to be counterbalanced by their lower signal-to-noise ratio. 

The two main sources of errors in the abundances are measurement error and analysis error 
caused by errors in the adopted model atmosphere parameters. The scatter of the abundances 
provided by individual lines of the same species, which are caused by measurement errors of the 
equivalent widths and, to a lesser extent, by errors in the line oscillator strengths, is a good guide 
to measurement error. This scatter, as measured by the standard deviation of the individual 
line abundances, is given in Table 3. (Although the standard deviations of the abundances from 
individual lines are larger than the standard deviations of the mean abundances, the heavy 
element abundances are based on only a few lines for each element — see Table 2 — so we prefer 
to consider the standard deviations of the individual line abundances.) For Rb, whose abundance 
is based on spectrum synthesis of the 7800 A Rb i line. Table 3 gives errors estimated from the 
fit of the observed and synthesised spectra. Inspection of Table 3 shows the measurement-related 
abundance errors range up to ±0.15 dex with larger errors in a few cases; a representative error in 
the [X/H] abundances is ~ ±0.07 dex, while a representative error of the [X/Fe] abundances is 
~ ±0.1 dex. 

Estimated errors in the adopted effective temperatures are between ±50 K, for stars with 
a good selection of Fei lines providing a well-determined excitation temperature, and ±100 K, 
for stars for which we adopted color-based effective temperatures. Representative errors in the 
adopted log g and metallicities are ±0.2 and ±0.1 dex, respectively. A representative uncertainty 
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in the microturbulence is ±0.5 km s , although we note that the Rb abundances provided by 
the weak Rb I hne have httle, or no, microturbulence dependence and that because of the metal 
deficiency of most of the program stars the abundances of most other elements also have only a 
small, or negligible, microturbulence dependence. Adopting a representative effective temperature 
error of ±100 K and the stated errors of the other parameters, we find that for a typical dwarf 
the combined effects of these errors change the Fe abundance (from Fei lines) by ±0.11 dex. The 
corresponding figure for a typical giant is ±0.17 dex. For the heavy elements the abundance of 
the element relative to Fe, [El/Fe], holds the most interest. This ratio is less dependent on the 
atmospheric parameters than the absolute abundance. In the typical dwarf the combined effects 
of the errors in the atmospheric parameters change [El/Fe] by ±0.04 (Rb), ±0.05 (Y), ±0.05 (Zr), 
±0.11 (Ba), and ±0.06 dex (Nd). The corresponding figures for a typical giant are: ±0.07 (Rb), 
±0.06 (Y), ±0.00 (Zr), ±0.14 (Ba), and ±0.07 dex (Nd). We estimate representative total errors, 
caused by measurement error and errors in the model atmosphere parameters together, to be 0.1 
- 0.2 dex in [Fe/H] and 0.1 - 0.2 dex for [El/Fe]. 

We now briefly consider how the Fe and Rb abundances of the dwarfs and subgiants would 
change if the preliminary effective temperatures, which are mostly based on the infrared flux 
method (Alonso et al. 1996a), and preliminary gravities, which are mostly Hipparcos-based, 
were used instead of the adopted effective temperatures and gravities. As discussed earlier, the 
preliminary effective temperatures were revised upward by an average of 140±75 K (cr of the 
individual differences) for 13 of the dwarfs and subgiants, while no revisions of the preliminary 
effective temperatures were made for the other 19 dwarfs and subgiants. Use of the lower 
preliminary effective temperatures for these 13 stars would decrease their Fe abundances (from 
Eel lines) by an average of 0.14dex, while their [Rb/Fe], as set by Rbi and Eel lines, would 
increase by an average of 0.04 dex. The adopted gravities for 28 of the dwarfs and subgiants are 
spectroscopically determined, while those of the other four dwarfs and subgiants, for which we 
could not determine spectroscopic gravities, are the preliminary gravities. Use of the preliminary 
gravities, instead of the adopted spectroscopic gravities, for these 28 stars would not change their 
Fe or Rb abundances significantly; the adoption of the preliminary log g, which are 0.12±0.12dex 
(a of the individual differences) higher on average than the spectroscopic log g, in place of the 
spectroscopic log g would change the [Fe/H] and [Rb/Fe] from neutral lines by ±0.01 and 0.00 dex, 
respectively, on average. 



4. RESULTS 

4.1. Comparison with the Literature 

Before we consider the Rb and other heavy-element abundances, we compare our results with 
those published in the literature. First we consider the results for [Fe/H]. Figure |2| compares the 
[Fe/H] determinations for stars in common to this study and earlier high signal-to-noise ratio. 



- 12 - 



high resolution studies. The comparison is not exhaustive, but does include all recent studies 
(since 1990) which have two, or more, stars in common with the present study. The agreement of 
the [Fe/H] determinations is good over most of the metallicity range and, although the [Fe/H] of 
this study tend to be slightly more negative than the literature [Fe/H] in the most metal-deficient 
stars, the overall agreement is not unsatisfactory. 

Previous studies which, to our knowledge, have determined Rb abundances for stars in 
common with those of the present investigation are Mackle et al.'s (1975a) study of Arcturus and 
Gratton & Sneden's (1994) study of heavy-element abundances in metal-poor stars. In Table 4 
we compare our [Rb/Fe] with those of the two earlier investigations. Because our Rb abundances 
are based on the 7800 A Rbi line, while Mackle et al.'s (1975a) abundance is based on both the 
7800 and 7947 A lines, we also include in the Table their abundance for the 7800 A Rbi line 
alone (Mackle et al. 1975b). (Gratton & Sneden's results are based only on the 7800 A line.) 
In order to make the Rb abundances of our and the earlier studies directly comparable we have 
adjusted the Rb abundances of the earlier studies to reflect the values they would have if the 
stellar parameters (Teg, log (7, [M/H], and ^) used in the earlier studies had been the same as 
those used here. For Arcturus the difference between our [Rb/Fe] and Mackle et al.'s is only -0.03 
dex — pleasingly small and not unexpected for the case of such a bright star. For the three stars 
that we have in common with Gratton & Sneden we note that both investigations determined 
Rb abundances for two of the stars (HD 64606 and 187111), but were only able to determine 
upper limits for the third star (HD 122956). The differences between our [Rb/Fe] and Gratton & 
Sneden's are -0.25 (HD 64606), -0.09 (HD 122956), and +0.17 (HD 187111). That these differences 
are much larger than in the case of Arcturus can be ascribed to the fact that all three stars 
are quite metal-deficient and have only a weak or undetectable Rb line. We estimate that for 
these three stars the uncertainty in the Rb abundance associated with fitting our observed and 
synthetic spectra of the Rbi line is 0.1 dex for HD 64606 and 187111 and 0.2 dex for HD 122956; 
Gratton & Sneden's Rb abundances must be subject to similar uncertainties also. We conclude, 
therefore, that our [Rb/Fe] and Gratton & Sneden's are probably the same to within the errors of 
measurement for these three stars. 



4.2. The Abundances of Rb, Y, Zr, Ba, and Nd 

As is customary, the abundances in Table 3 are plotted as [el/Fe] against [Fe/H], in Fig. ^, 
^, and |5| to reveal trends in the relative abundance of element el and iron. Two points are 
immediately apparent: (i) three stars are unusually rich in the heavy elements — note especially 
the [Y/Fe] ratios of HD 23439A and B, and BD -^5° 3640 which we shall dub CH stars, and (ii) 
the distinctive behavior of [Rb/Fe] in the metal-poor stars — [Rb/Fe] > when the other heavy 
elements show [el/Fe] ~ 0. Before commenting on these striking results, we compare our results 
for the Y, Zr, Ba, and Nd abundances with results in the literature. 

Previous extensive abundance determinations of heavy (and other) elements in metal-poor 
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stars have shown that the run of [el/Fe] against [Fe/H] is smooth down to about [Fe/H] = -2 
with a 'cosmic' scatter less than the scatter that results from the errors of measurement. Cosmic 
scatter is present for more metal-poor stars but our sample is devoid of such stars. Therefore, our 
results are expected to agree well with previous studies despite the lack of a complete overlap in 
stellar samples. Key papers reporting results on heavy elements are Zhao &: Magain (1991) and 
Gratton & Sneden (1994) with reviews by Wheeler, Sneden, & Truran (1989), Lambert (1989), and 
McWilliam (1997) amongst others. Our results for [Y/Fe], [Ba/Fe], and [Nd/Fe] are in excellent 
agreement with previous results, for example, Zhao & Magain (1991) and Gratton & Sneden 
(1994) find [Y/Fe] ~ —0.1 at [Fe/H] = -1 with the relative underabundance of Y increasing to 
about 0.25 at [Fe/H] = -2 which agree well with Figure]^. A discrepancy appears when comparing 
results for [Zr/Fe]. Zhao & Magain (1991) and Gratton & Sneden (1994) report [Zr/Fe] ~ +0.2 
for [Fe/H] in the range of -1 to -2 but our results (Figure Q) show [Zr/Fe] to be consistently less 
than zero: a difference in [Zr/Fe] of about 0.3 to 0.4 dex relative to the previous studies. This 
difference is most probably due to our exclusive use of Zri lines. Brown, Tomkin, & Lambert 
(1983) found that Zri lines in mildly metal-poor giants gave a clear Zr underabundance which 
was plausibly attributed to non-LTE effects such as over-ionization of Zr atoms to Zr"*" ions. For 
metal-poor stars, Gratton & Sneden remark that their selection of Zr i lines gives a systematically 
lower Zr abundance than the Zr ii lines: the difference of -0. 16+0.05 dex would account in part for 
our largely negative values of [Zr/Fe]. 

Our results clearly show a relative overabundance of Rb in metal-poor stars: the mean value 
[Rb/Fe] = +0.23+0.02 (s.e.) is found from nine stars with [Fe/H] < —1, excluding the three CH 
stars. This is consistent with the four measurements reported by Gratton &; Sneden (1994) and 
with their six upper limits to [Rb/Fe]. 

Non-LTE effects such as overionization warrant consideration. As a guide to the non-LTE 
effects on Rb, we consider those calculated for lithium, another alkali. Carlsson et al. (1994) 
predict that the LTE abundances from the Lii 6707 A resonance doublet require correction by 
not more than 0.03 dex for non-LTE effects, a negligible correction in the present and almost all 
contexts. The abundances of Li and Rb are similar: lithium has the abundance log e(Li) ~ 2.2 
in the warmer dwarf stars comprising the Spite plateau and Rb declines from log e(Rb) = 2.6 at 
solar metallicity to 1.9 at [Fe/H] = -1 and to 0.8 at [Fe/H] = -2. The key point is that optical 
depth effects in lines and continua are slight for both elements. Rb is probably more affected by 
photoionization because the atom's ionization potential is 4.18 eV versus 5.39 eV for the lithium 
atom. Photoionization of Rb will be enhanced relative to the rate for Li but collisional ionization 
rates will also be enhanced. The different wavelengths of the resonance (and excited) atomic 
lines for Rb and Li will introduce no more than slight differences in the non-LTE corrections. 
If non-LTE effects were large for Rb, we would anticipate that dwarfs and giants of the same 
metallicity would yield systematically different Rb abundances. This is not the case: four dwarfs 
with [Fe/H] in the range -1 to -2 give a mean [Rb/Fe] = +0.28+0.04 (s.e.) and five giants in the 
same [Fe/H] range give the mean [Rb/Fe] = +0.18+0.02 (s.e.). (The mean for the dwarfs excludes 
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the three CH stars.) Although a non-LTE analysis for Rb would be of interest, we suggest that 
our results derived from LTE analyses are not substantially different from non-LTE results. 

Scatter of the [el/Fe] results at a given [Fe/H] is not significantly different from that expected 
from the measurement errors. Obviously, the three stars over-abundant in the heavy elements and 
dubbed CH stars are set aside as special cases. The scatter for [Rb/Fe] between [Fe/H] of -0.5 and 
-1.0 is small and consistent with the measurement errors. There is an apparent moderate increase 
in scatter of [Rb/Fe] below [Fe/H] = -1 but this is probably again due to the measurement errors 
because the Rbi line is very weak in these metal-poor stars. The results are roughly consistent 
with a constant [Rb/Fe] in stars with [Fe/H] < — 1. 

4.3. The New CH Stars HD 23439A and B and BD -^5° 3640 

These three stars, which are all dwarfs, are consistently overabundant in all of the five heavy 
elements investigated in this study. Figure ^ shows a Zr i line and a V i line in HD 23439A and B 
and HD 103095, a non-CH star with otherwise similar properties. The much greater strength of 
the Zr I line relative to the Vi line in HD 23439A and B as compared with HD 103095 is evident. 
As may be seen in Fig. ^, ^, and |^ the heavy-element abundance enhancements are similar in the 
three stars. The average enhancements for the three stars are: [Rb/Fe] = -1-0.41, [Y/Fe] = -1-0.34, 
[Zr/Fe] = +0.51, [Ba/Fe] = -^0.27, and [Nd/Fe] = -^0.32 (HD23439A and BD-F5°3640); actual 
values of [X/Fe] for the individual stars are given in Table 5. HD 23439 is a nearby visual binary 
composed of a KIV primary and a K2V secondary; the Hipparcos Catalogue gives 7."307 and 
40.83 mas for the separation of its components and parallax, respectively. These numbers set a 
lower limit on the A-B linear separation of 179 AU. To the best of our knowledge HD 23439 is the 
first case of a binary in which both components have been found to be CH stars. 

We note that HD 23439B is a single-lined spectroscopic binary with a period of 48.7 d and 
a mass function of 0.0022 (Latham et al. 1988). Could the unseen companion of HD23439B 
be the white dwarf descendant of the AGB star responsible for the mass transfer that changed 
HD 23439 A and B into CH stars? Perhaps, but in this scheme it is hard to explain the very 
similar heavy-element enhancements of HD23439A and B (see Table 5). Just as today the unseen 
companion is much closer to component B than to component A so in the past the putative AGB 
predecessor of the unseen companion must also have been much closer to B than A. How did 
the AGB star manage to give the A and B components the same heavy-element enhancements? 
This difficulty with the AGB-star scenario suggests that the heavy-element enhancements must be 
primordial. 
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5. RUBIDIUM AND STELLAR NUCLEOSYNTHESIS 

Heavy elements are synthesised by the neutron capture s- and r-processes. (In considering 
elemental abundances, the small contribution from p-processes may be neglected.) Detailed 
dissection of the isotopic abundances measured for carbonaceous chondrites has provided an 
isotope by isotope resolution of the abundances into s- and r-process contributions (cf. Kappeler, 
Beer, & Wisshak 1989). As is well known, Ba and Eu are primarily s- and r-process products 
respectively: Cowan (1998) estimates that Ba is 85% an s-process product, and Eu is 97% a 
r-process product — see Gratton & Sneden (1994) for similar estimates. Rubidium is of mixed 
parentage; Cowan gives the s- and r-process fractions as 50% each — Gratton & Sneden provide 
quite similar estimates (48% for s- and 52% for the r-process). As noted in the Introduction, 
the s-process contribution may be broken into a 'weak' and a 'main' component. Gratton & 
Sneden divide the 48% s-process Rb contribution into 5% from the weak and 43% from the main 
s-process. Our goal is to use the Ba and Eu abundances as monitors of the s- and r-processes 
respectively to predict the Rb abundances, and then to comment on the consistency between the 
predicted and observed Rb abundances. 

The other heavy elements considered here are also a mix of s- and r-processes: Cowan gives 
the following (s,r) %: (72, 28) for Y, (81, 19) for Zr, and (47, 53) for Nd. Gratton & Sneden 
(1994) put the weak contribution to the total s-process as 16% for Y, 10% for Zr, and less than 
1% for Ba, Nd, and Eu. Especially interesting is the roughly 50-50 split for Nd that matches the 
split for Rb. Then, the simplest possible scenario of unvarying yields of s and r-process products 
over the life of the Galaxy would predict that [Rb/Fe] and [Nd/Fe] would vary identically with 
[Fe/H]. Inspection of Fig. ^ and |5| shows that this is not the case. 

There are several factors pertinent to the understanding of the run of Rb and other heavy 
elements with [Fe/H]. 

• It is now well known that, the distribution of heavy elements at low [Fe/H] resembles a 
r-process pattern and is not the mix of s- and r-processes that prevails at solar metallicities 
(Truran 1981; Sneden & Parthasarathy 1983; Wheeler et al. 1989; Lambert 1989). There 
is evidence that the abundance distribution of the r-process was largely invariant from low 
metallicities to the present [Fe/H] ~ 0. Then, it should suffice in modelling the [el/Fe] vs 
[Fe/H] relations to adopt the relative r-process abundances that are obtained from dissection 
of the measurements on carbonaceous chondrites. 

• Europium is assigned to the r-process: Cowan's (1998) resolution of the meteoritic 
abundances is 97% r-process and a mere 3% s-process. With declining metallicity the s/r 
ratio declines. Therefore, we may assume that Eu is a r-process product throughout the 
evolution of the Galaxy. The run of [Eu/Fe] against [Fe/H] is taken from McWilliam's (1997) 
review: [Eu/Fe] = at [Fe/H] = with a smooth transition to [Eu/Fe] ~ 0.3 at [Fe/H] = 
-1.0 and to the metallicity limit [Fe/H] = -2.5 of interest to us. 
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• Red giants enriched in s-process heavy elements are Hkely the major donors of these elements 
to the Galaxy's interstellar medium. Analyses of such red giants of differing [Fe/H] show 
that the pattern of s-process products has evolved with [Fe/H]. Smith (1997), who has 
collated published results, defines Y and Zr as 'light' s-process elements (here, Is) and Ba, 
La, and Ce as 'heavy' s-process elements (here, hs). He finds that [hs/ls], which by definition 
is at [Fe/H] = 0, increases to [hs/ls] ~ 0.6 at [Fe/H] = -1.5. This evolution of [hs/ls] is 
attributed to an increase in the average exposure to neutrons in the He-burning shell of the 
AGB stars that are the site of the main s-process, as expected on theoretical grounds. 

• There is limited evidence also from abundance analyses of red giants that the Rb abundance 
relative to other Is elements increases with decreasing [Fe/H]; Smith's (1997) collection of 
results implies [Rb/Zr] increases by about 0.7 dex from [Fe/H] = to -1.5. This increase is 
attributed to a higher mean neutron density in the He-burning shell of the metal-poor AGB 
stars. As noted by Smith, this increase is expected on theoretical grounds. 

• If the magnitude of the weak s-process contributions to the abundances in carbonaceous 
chondrites were representative of the s-process at all relevant [Fe/H], the weak s-process 
could be safely dropped from our search for an explanation of the run of [el/Fe]. A thorough 
direct check on the weak s-process is not possible because the majority of the elements 
between Rb and the Fe-group are inaccessible spectroscopically. Zinc, which offers at 
least a hint of the behavior of the weak s-process, is assigned 34% to the s-process and 
66% to the r-process in Cowan's breakdown of solar system abundances. The s-process 
component is essentially entirely due to the weak s-process. This breakdown neglects a 
possible contribution to Zn from the sources that contribute the Fe-group elements. In the 
metallicity range of interest, [Zn/Fe] = 0.0ib0.15 (Sneden & Crocker 1988). This result, 
which is barely compatible with the increase of [Eu/Fe] with declining metallicity, implies a 
drop in [s/Fe] with metallicity and justifies our neglect of the weak s-process contribution to 
Rb and other elements. 

Guided by these facts, it is possible to predict relative abundances of the heavy elements 
including Rb. We begin by considering Ba, Nd, and Eu. Eu defines the evolution of the r-process 
products. The observed run of [Ba/Fe] against [Fe/H] provides the evolution of the hs component 
of the s-process after a small correction for this element's r-process component based on the Eu 
abundances and the meteoritic r-process Ba/Eu ratio. The adopted runs of [Eu/Fe] and [Ba/Fe] 
against [Fe/H] are shown in Fig. Then, it is a simple matter to predict the run of [Nd/Fe] using 
the meteoritic 50-50 split into s- and r-process contributions. This prediction which is shown 
too in Fig. ^a is slightly inconsistent with the observations (Fig. H) that show [Nd/Fe] ~ at all 
metallicities rather than the predicted [Nd/Fe] = 0.14 at [Fe/H] < ~1. Earlier, we noted that our 
Nd abundances are quite consistent with previously published results. This small inconsistency 
appears not to have been noted previously. It is likely that when the measurement errors are 
included the prediction and observations will overlap. Note that the prediction uses the Eu and 
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Ba abundances as well as the meteoritic s to r ratios for Ba, Nd, and Eu. A change in the Nd 
ratio from 50% s and 50% r to 75% s and 25% r reduces [Nd/Fe] to 0.0 for metal-poor stars. 

Yttrium abundances may be predicted from the [Ba/Fe] observations and Smith's estimates 
of [hs/ls] from heavy element enriched red giants such as S and Barium stars. This prediction 
is shown in Fig. |^. The relative underabundance of Y (i.e., [Y/Fe] < 0) results largely from 
the steep increase in [hs/ls] with decreasing metallicity that offsets the increase in the r-process 
contribution. This particular prediction assumes a meteoritic ratio of 72% s-process and 28% 
r-process and makes no attempt to separate main from weak s-process contributions. This 
prediction matches the observations quite well (see Fig. ^). 

The Rb prediction corresponding to the Y prediction which is also shown in Figure 03 does 
not correspond to the observed [Rb/Fe] ratios in metal-poor stars. The limited evidence gathered 
by Smith suggests that Rb in heavy element enriched red giants is progressively overabundant 
(relative to Zr) in metal-poor giants. This increase is attributable to a higher neutron density 
at the s-process site in red giants. If a smooth curve is drawn through Smith's assembled data, 
the resulting run of [Rb/Fe] is shown in Fig. |^. This prediction is in good agreement with the 
observations. 

6. CONCLUDING REMARKS 

The principal novel result of our survey of Rb abundances in stars is that rubidium relative to 
iron is systematically overabundant in metal-poor stars. This increase reflects partly the growth 
of the r-process abundances relative to iron in metal-poor stars; we model this increase using 
observed Eu abundances and the assumption that the pattern of r-process abundances is solar-like 
at all metallicities. A second and major factor accounting for the increase in the Rb to Fe ratio in 
metal-poor stars is that the s-process contribution to Rb increases with decreasing metallicity. We 
model this increase using published abundances of Rb and other heavy elements collated by Smith 
(1997) for s-process enriched red giants that are presumed to be representative of the donors of 
s-processed material to the interstellar medium and so to control the chemical evolution of the 
Galaxy as far as the s-process is concerned. Two factors influence the Rb abundance: (i) the total 
exposure to neutrons at the s-process site, the He-burning shell of an AGB star, increases with 
decreasing metallicity of the red giant, and (ii) the neutron density at the s-process site increases 
with decreasing metallicity. That (i) is true follows from the observed increase of the relative 
abundance ratio of heavy to light s-process elements in the s-process enriched red giants. It is 
this effect that accounts, for example, for the drop in [Y/Fe] in metal-poor stars. That (ii) is true 
follows from the limited data on Rb abundances in s-process enriched red giants. As explained 
above, neutron density influences the Rb abundance through the branch in the s-process path at 
^^Kr. In this thoroughly empirical way we account for the relative enrichment of Rb in metal-poor 
stars. In short, the observed [Rb/Fe] ratios of metal-poor stars are consistent with the expectation 
that AGB stars control the input of main s-process products to the Galaxy's interstellar medium. 
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A serendipitous discovery is the finding that both members of a visual binary are mild 
CH stars, HD23439A and B, with s-process overabundances relative to other stars of the same 
metallieity. Mass transfer across a binary system is now thought to account for CH stars and the 
Barium stars, the higher metallieity counterparts of the CH stars. HD 23439B is a single-lined 
spectroscopic binary and the visible star might have been transformed to a CH star by mass 
transfer from the companion, then an AGB star and now a white dwarf. HD 23439A appears to 
be a single star that cannot have captured significant amounts of mass from a very distant AGB 
star orbiting HD 23439B. We suggested, therefore, that these CH stars testify that the halo's 
interstellar medium was not entirely chemically homogeneous. This is not too surprising given 
that s-process products are injected into the interstellar medium at low velocity by red giants 
whereas iron and other elements are injected at very high velocity by supernovae. If the timescale 
for star formation is shorter than the timescale for thorough mixing of supernovae and red giant 
ejecta, abundance anomalies will result. 

We thank Bengt Edwardsson for providing the grid of dwarf and subgiant MARCS model 
atmospheres, Chris Sneden for providing Pilachowski et al.'s (1996) giant MARCS model 
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research has made use of the Simbad database, operated at CDS, Strasbourg, Prance. This work 
has been supported in part by NSF grant AST 9618414 and the Robert A. Welch Foundation of 
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Figure Captions 

Fig. 1. — The Rb I line at 7800.27 A and the adjacent Si I hne at 7800.00 A in three of the program 
stars. The top panel shows a 200,000 resolving-power observation of jj, Cas (G5Vb) and the middle 
and bottom panels show 60,000 resolving-power observations of HD 65583 (G8V) and HD 108564 
(K2V), respectively. Synthesised spectra for [Rb/Fe] of 0.0, 0.3, and 0.6 are shown. In ji Cas and 
HD 65583 the high-excitation Sil line (xiower = 6.18 eV) is much stronger than the low-excitation 
Rbl line (xiowcr = 0.00 eV), but in HD 108564, which is ~700K cooler than Cas and HD 65583, 
the Rb I line is slightly stronger than the Si I line. 

Fig. 2. — Comparison of derived [Fe/H] for stars in common to the present study and earlier high 

signal-to-noise ratio, high resolution studies. Literature shorthand citations are: P96 = Pilachowski 
et al. (1996); S91 = Snedcn, Gratton, & Crocker (1991); F98 = Fuhrmann (1998); F97 = Flynn 
& Morell (1997); D93 = Drake & Smith (1993); T92 = Tomkin et al. (1992); P93 = Pilachowski, 
Sneden, & Booth (1993); G94 = Gratton & Sneden (1994); K97 = King (1997); G99 = Gratton et 
al. (1999). The dashed line shows the locus of identical results. 

Fig. 3. — Plots of [Rb/Fe] and [Y/Fe] versus [Fe/H] for the program stars. The three CH stars 
(HD23439A and B and BD-K5°3640) are shown with filled-in squares. The error bars for [Rb/Fe] 
are a quadratic sum of the Rb abundance errors, estimated from the fit of the observed and synthetic 
spectra of the 7800 A Rb l line, and the standard deviations of the Fe abundances from individual 
lines. The error bars for [Y/Fe] are a quadratic sum of the standard deviations of the Y and Fe 
abundances from individual lines. 

Fig. 4. — Plots of [Zr/Fe] and [Ba/Fe] versus [Fe/H] for the program stars. The three CH stars 
(HD23439A and B and BD -|-5° 3640) are shown with filled-in squares. The error bars for [Zr/Fe] 
are the quadratic sum of the standard deviations of the Zr and Fe abundances from individual lines. 
Likewise for [Ba/Fc]. 

Fig. 5.— Plot of [Nd/Fe] versus [Fe/H] for the program stars. The two CH stars (HD23439A 
and BD +5° 3640) with Nd abundance determinations are shown with filled-in squares. The Nd 
abundances are based on a single Ndii line so a fixed error bar is shown for [Nd/Fe]. 

Fig. 6. — A Zrl line and a Vl line in the two CH stars HD 23439 A and B and in a comparison star 
(HD 103095). The Zr l and V l lines have similar excitation potentials, 0.00 and 1.05 eV, respectively, 
so their relative strengths are set primarily by the relative Zr and V abundances. The much greater 
strength of the Zr I line relative to the V I hne in HD 23439A and B compared with HD 103095 
reveals the enhancement of Zr in HD23439A and B. 

Fig. 7. — a) The runs of the "heavy" neutron-capture elements Eu, Nd, and Ba as a function of 
[Fe/H]. The [el/Fe] for Eu and Ba, which are primarily r- and s-process elements, respectively, are 
schematic representations of the observed [el/Fe], while that for Nd is predicted on the basis of the 
behavior of Eu and Ba. See text for details, b) The predicted runs of the "light" neutron-capture 
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elements Rb and Y as a function of [Fe/H]. These predicted abundances for Rb and Y, which also 
rest on the behavior of Eu and Ba, include allowance for the evolution of [hs/ls] with [Fe/H]. Two 
predictions are shown for Rb. See text for details. 
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Table 1. Rb Program Stars, Model Atmosphere Parameters 



Name/HD/DM 


V 




Sp. T. 


[Fe/H] 




log 


g 


[M/H] 




Dwarfs and Subgiants 


// Cas 


5, 


.16 


G5Vb 


-0.92 


5350 


4 


.47 


-0 


.91 


1.0 


r Get 


3 


.50 


G8V 


-0.59 


5330 


4 


.30 


-0 


.57 


0.5 


61 CygA 


5. 


.22 


K5V 


-0.43 


4450 


4 


.56 


-0 


.33 


0.5 


61 CygB 


6 


.03 


K7V 


-0.63 


4120 


4 


.40 


-0 


.50 


1.0 


e Eri 


3 


.73 


K2V 


-0.14 


5110 


4 


.38 


-0 


.14 


0.7 


4906 


8, 


.77 


GO 


-0.75 


5150 


3 


.52 


-0 


.74 


1.0 


18907 


5, 


.88 


G5IV 


-0.74 


5030 


3 


.38 


-0 


.72 


1.0 


23439A 


8. 


.19 


KIV 


-1.05 


5200 


4 


.61 


-1 


.05 


0.5 


23439B 


8. 


.78 


K2V 


-1.11 


4950 


4 


.71 


-1 


.10 


1.0 


24341 


7, 


.88 


GIV 


-0.69 


5450 


3 


.59 


-0 


.68 


1.0 


25329 


8, 


.49 


KIV 


-1.82 


4775 


4 


.71 


-1 


.81 


1.0 


45282 


8, 


.03 


GO 


-1.58 


5275 


3 


.00 


-1 


.57 


1.0 


64606 


7. 


.43 


G8V 


-0.95 


5245 


4 


.35 


-0 


.93 


1.0 


65583 


6 


.97 


G8V 


-0.82 


5305 


4 


.44 


-0 


.80 


1.0 


103095 


6 


.42 


G8Vp 


-1.48 


5040 


4 


.31 


-1 


.44 


1.0 


108564 


9 


.43 


K2V 


-1.18 


4660 


4 


.67 


-1 


.06 


1.0 


134439 


9 


.06 


KO 


-1.54 


5040 


4 


.26 


-1 


.48 


1.0 


134440 


9 


,42 


KOV: 


-1.44 


4945 


4 


.72 


-1 


.42 


1.0 


148408 


9 


.64 


G5 


-0.84 


5355 


4 


.-12 


-0 


.83 


1.0 


149414 


9 


.61 


G5V 


-1.36 


5180 


4 


.27 


-1 


.35 


0.5 


150281 


8, 


.65 


KO 


-0.02 


5150 


4 


.26 


-0 


.07 


0.5 


192031 


8, 


.67 


G8V 


-0.89 


5350 


4 


.54 


-0 


.88 


1.0 


212753 


9 


.67 


G5 


-0.51 


5100 


4 


.50 


-0 


.48 


0.5 


216179 


9 


.33 


G 


-0.86 


5395 


4 


.38 


-0 


.86 


1.0 


224618 


8, 


.92 


KOV 


-0.76 


5320 


4 


.50 


-0 


.74 


1.0 


230409 


10 


.10 


GO 


-1.00 


5320 


4 


.-10 


-0 


.98 


1.0 


233832 


10 


.14 


KOV 


-0.83 


4970 


4 


.37 


-0 


.78 


1.0 


-22°3557 


10 


.27 


KO 


-0.87 


5055 


4 


.60 


-0 


.83 


1.0 


+5°3640 


10 


.35 




-1.27 


5080 


4 


.39 


-1 


.26 


1.0 


+24<'4460 


9 


.47 


KOV 


-0.89 


5185 


4 


.38 


-0 


.84 


1.0 


+41 "3306 


8 


.87 


KOV 


-0.62 


5160 


4 


.-13 


-0 


.58 


1.0 


+48°3755 


9 


.57 


G8 


-0.63 


5325 


4 


.47 


-0 


.63 


1.0 


Giants 


Arcturus 


-0.04 


K1.5III 


-0.63 


4300 


1.50 


-0.50 


1.7 


3008 


9.7 


KO 


-1.95 


4300 


0.60 


-2.00 


2.2 


6833 


6.75 


G9III 


-0.80 


4580 


1.86 


-0.80 


1.7 


8724 


8.30 


G5 


-1.92 


4500 


1.02 


-1.91 


1.6 


29574 


8.33 


G8/K0III 


-1.95 


4280 


0.22 


-1.80 


1.7 


74462 


8.74 


G5IV 


-1.61 


4600 


1.68 


-1.61 


1.8 


83212 


8.28 


G8III 


-1.57 


4550 


1 


.14 


-1 


.59 


1.8 


122956 


7.22 


G6IV 


-1.88 


4600 


1.20 


-1 


.91 


1.8 


187111 


7.72 


G8 


-1.72 


4400 


1.08 


-1.70 


1.7 


204543 


8.6 


GO 


-1.89 


4700 


1.28 


-1.93 


2.2 


220838 


9.30 


G8III 


-1.94 


4300 


0.84 


-1.97 


2.2 


-9°5831 


9.5 




-1.94 


4575 


1.12 


-1.80 


1.5 
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Table 2. Atomic Line Data and Solar Equivalent Widths 



A 


\ 


logs/ 


Source of gf 




Waq 


(A) 


(eV) 








(mA) 


FeT logfg = 7.51 








5849.69 


3.69 


-2.84 


MRW 


1.4 


7.8 


5852.23 


4.55 


-1.18 


MRW 


1.4 


41.2 




4.D1 


1 /I Q 
— 1 .4o 


RKK 


1.4 


22.4 


5856.10 


4.29 


-1.57 


solnr 


1.4 


33.4 


5858.79 


4.22 


-2.11 


MRW 


1.4 


13.8 


5859.60 


4.55 


-0.59 


solar 


1.4 


72.8 


5861.11 


4.28 


-2.33 


solar 


1.4 


8.9 


5862.37 


4.55 


-0.35 


solar 


1.4 


fSfS.4 


6151.62 


2.18 


-3.31 


O'B, B, BKK 


1.4 


50.9 


6157.73 


4.08 


-1.11 


MRW 


1.4 


63.4 


6159.38 


4.61 


-1.86 


solar 


1.4 


12.6 


6165.36 


4.14 


-1.47 


O'B 


1.4 




6173.34 


2.22 


-2.88 


B 


1.4 


69.7 


6591 .33 


4.59 


-1.97 


solar 


1.4 


10.8 


6608 04 


2.28 


-3.88 


MRW 


1.2 


18.7 


6713.75 


4.80 


-1.45 


MRW 


1.4 


21.5 


6725.36 


4.10 


-2.17 


solar 


1.4 


17.9 


6733.15 


4.64 


-1.43 


MRW 


1.4 


27.7 


6857.25 


4.08 


-2.00 


MRW 


1.4 


23.5 


7802 51 


5.09 


-1.34 


solar 


1.4 


16.1 


Fe II loj 




7.51 








5425.26 


3.20 


-3.29 


solar 


2.5 


41 .8 


6149.25 


3.89 


-2.76 


solar 


2.5 


37.0 


6247.56 


3.89 


-2.38 


solar 


2.5 


54.4 


6369.46 


2.89 


-4.15 


solar 


2.5 


19.9 


6456.39 


3.90 


-2.17 


solar 


2.5 


64.7 


Rb I lof 


5^0 = 


2.60 








7800.29 


0.00 


+0.13 


WM 


2.5 




Yn log 


fg = 2.24 








4883.69 


1.08 


+0.07 


H 


2.5 


60.0 


4900.12 


1.03 


-0.09 


H 


2.5 


57.9 


5087.43 


1.08 


-0.17 


H 


2.5 


48.3 


Zr I log 


£q = 2.60 








6127.48 


0.15 


-0.60 


solar 


2.5 


3.8 


6134.57 


0.00 


-0.89 


solar 


2.5 


2.8 


6143.18 


0.07 


-0.72 


solar 


2.5 


3.5 


Ball lo 


8^0 = 


2.13 








5853.69 


0.60 


-1.01 


WM 


3.0 


64.2 


6141.73 


0.70 


-0.08 


WM 


3.0 


119.5 


6496.91 


0.60 


-0.38 


WM 


3.0 


102.7 
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Table 2. Atomic Line Data and Solar Equivalent Widths, continued 



X 


X logs/ 


Source of gf 




Waq 


(A) 


(eV) 






(mA) 


Ndll 


log fQ = 1.50 








5319.82 


0.55 -0.21 


W 


2.5 


11.4 



Notes lo Table 2. 

Solar abinidaTices are adopted from Anders & Grevesse 
(1989), except the Fe abundance which is from Lambert el: al. 
(1996). The columns give: 1) the line rest wavelength, 2) the 
excitation energy of the lower level of the transition, 3) the 
logarithm of the gf-value, 4) the source of the 5/- value. 5) an 
enhancement factor to the classical van der Waals damping 
constant. 6) the solar equivalent width from light reflected 
by the asteriod Iris. Sources of ^i/- values: B — Blackwell et 
al. (1976), see Fuhr, Martin, k Wiese (1988) for additional 
references; BKK — Bard, Kock, k Kock (1991); H — Han- 
naford et al. (1982); MRW — May, Richter, k Wichelmann 
(1974), increased by 0.11 dex if log ^i/ > —0.75 and 0.05 dex 
if log 3/ < -0.75; O'B — O'Brian et al. (1991); solar — 
solar ^i/- values calculated with the solar abundances and so- 
lar equivalent widths given and Holweger k Muller's (1974) 
model atmosphere; W — Ward et al. (1985); WM — Wiese 
k Martin (1980). 
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TABLE 3 

Abundance Results, [X/H] 



Name/HD/DM Fel Fell Rbl Yn Zrl Ball Ndll 



Dwarfs and Subgiants 
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n 09 

— U.yZ 






— U.UU zt 


1 
. lU 






— 1 .UD 












65583 


-0.82 ± 


.06 


-0.80 






-0.62 ± 


.10 






-0.96 ± 


.03 










103095 


-1.48 ± 


.03 


-1.46 ± 


.05 




-1.40 ± 


.20 


-1.65± 


.05 


-1.46 ± 


.02 


-1.55± 


.03 






108564 


-1.18± 


.13 








-0.98 ± 


.08 






-1.22 ± 


.10 










134439 


-1.54 ± 


.03 


-1.52 




< 


-1.34 ± 


.20 


-1.99 ± 


.10 






-1.77± 


.04 






134440 


-1.44 ± 


.08 


-1.45 




< 


-1.24± 


.30 


-1.81± 


.09 






-1.62± 


.04 






148408 


-0.84 ± 


.06 


-0.85 ± 


.09 




-0.49 ± 


.10 


-0.93 ± 


.09 






-0.93 ± 


.08 


-0 


77 


149414 


-1.36± 


.05 


-1.33 ± 


.04 




-0.91 ± 


.15 


-1.21 








-1.24± 


.09 






150281 


-0.02 ± 


.06 


-0.02 ± 


.03 




-0.12 ± 


.08 


-0.16± 


.08 


-0.25 ± 


.05 


0.00 ± .09 


-0 


11 


192031 


-0.89 ± 


.06 


-0.89 ± 


.00 




-0.54 ± 


.10 


-0.88 ± 


.05 


-0.73 ± 


.05 


-0.85 ± 


.09 






212753 


-0.51 ± 


.09 


-0.50 ± 


.08 




-0.71 ± 


.10 


-0.68 ± 


.00 


-0.83 ± 


.22 


-0.41 ± 


.03 


-0 


42 


216179 


-0.86 ± 


.05 


-0.83 ± 


.03 




-0.56 ± 


.10 


-0.97 ± 


.04 


-0.91 ± 


.20 


-0.97± 


.05 


-0 


81 


224618 


-0.76 ± 


.06 


-0.76 ± 


.05 




-0.69 ± 


.10 


-0.88 ± 


.03 


-0.74 ± 


.03 


-0.79 ± 


.02 






230409 


-1.00 ± 


.06 


-0.96 ± 


.03 


< 


-1.00 ± 


.30 


-1.12± 


.04 






-1.08± 


.05 






233832 


-0.83 ± 


.04 


-0.81 






-0.58 ± 


.10 






-0.90 ± 


.06 


-0.96 








-22°3557 


-0.87 ± 


.06 


-0.85 ± 


.03 




-0.77 ± 


.10 


-l.OOit 


.06 


-0.85 ± 


.07 


-l.Oli 


.06 






+5°3640 


-1.27± 


.05 


-1.23 ± 


.02 




-0.74 ± 


.15 


-0.85 ± 


.07 


-0.68 ± 


.07 


-0.85± 


.09 


-0 


84 


+24M460 


-0.89 ± 


.07 


-0.86 ± 


.07 




-0.63 ± 


.10 


-0.95 ± 


.04 


-0.87 ± 


.09 


-0.95 ± 


.03 


-0 


96 


+41°3306 


-0.62 ± 


.07 


-0.61 ± 


.07 




-0.45 ± 


.10 


-0.79 ± 


.09 


-0.74 ± 


.05 


-0.77 ± 


.11 






+48°3755 


-0.63 ± 


.07 


-0.66 ± 


.03 




-0.35 ± 


.10 


-0.84 ± 


.03 


-0.58 




-0.68 ± 


.05 


-0 


52 


Giants 


Arcturus 


-0.63 ± 


.09 


-0.66 ± 


.04 




-0.58 ± 


.02 


-1.21 




-1.00 ± 


.01 


-0.67± 


.01 


-0 


63 


3008 


-1.95 ± 


.05 


-1.96 ± 


.01 




-1.75 ± 


.20 


-2.26 




-2.26 ± 


.04 


-2.05 ± 


.01 


-1 


95 


6833 


-0.80 ± 


.03 


-0.83 ± 


.04 




-0.84 ± 


.05 






-1.02± 


.01 


-0.81 








8724 


-1.92± 


.16 


-1.92 ± 


.08 




-1.69± 


.20 






-2.26 ± 


.21 


-1.85 








29574 


-1.95 ± 


.08 


-1.96 






-1.63 ± 


.20 






-2.11 ± 


.07 










74462 


-1.61± 


.10 


-1.62 ± 


.10 


< 


-1.41 ± 


.20 










-1.77 








83212 


-1.57± 


.15 


-1.56 






-1.47 ± 


.20 






-1.61 ± 


.04 










122956 


-1.88 ± 


.11 


-1.86 




< 


-1.48 ± 


.20 


















187111 


-1.72± 


.05 


-1.73 ± 


.04 




-1.67± 


.10 


-1.84± 


.09 


-1.91 ± 


.05 


-1.64± 


.01 


-1 


72 


204543 


-1.89 ± 


.08 


-1.84± 


.02 


< 


-1.49 ± 


.20 


-2.18± 


.06 






-1.68± 


.11 


-1 


88 


220838 


-1.94± 


.03 


-1.92 ± 


.06 


< 


-1.84± 


.25 


-2.24 




-2.05 ± 


.08 


-2.15± 


.01 


-1 


93 


-9°5831 


-1.94± 


.04 


-1.93 ± 


.02 


< 


-1.24± 


.30 


-2.26 ± 


.04 






-2.05 ± 


.09 


-2 


09 



Note. — The abinidanees of all elements, with the exception of Tib, are based on eqnivalent widths. The errors are 
the standard devialions of iho individual line abundances of the species; absence of an error means only one line 
was measured. The Tib abundances are based on spectrum synthesis of the 7800.3 A Rbl line. Errors for Rb are 
estimated from the fit of the observed and synthesised spectra. 
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TABLE 4 

Comparison of Rb Abundances with Literature 



Star 


[Fe/H] 




[Rh/Fe] 






Difference 


Reference and Lines Used 






this 


lit. 


adj. 




this study 








study 




lit. 




— adj. lit. 




Arcturus 


-0.63 


+0.05 


-0.05 


-0 


06 


+0.11 


Mackle et al. (1975a), 7800 and 7947 A 


Arcturus 


-0.63 


+0.05 


+0.09 


+0 


08 


-0.03 


Mackle et al. (1975b), 7800 A 


Hn 64606 


-O.flS 


+0.35 


+0.60 


+0 


60 


-0.25 


Grnlion k Sncdcn (1994). 7800 A 


IID 122956 


-1.88 


< +0.40 


< +0.40 


< +0 


49 


-0.09 


Gratton k Sneden (1994), 7800 A 


HD 187111 


-1.72 


+0.05 


0.00 


-0 


12 


+0.17 


Gratton k Sneden (1994), 7800 A 



1 
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TABLE 5 

Heavy-Element Abundances oe the CH Stars 



HD/DM 


[Fe/H] 


[Rb/Fe] 


[Y/Fe] 


[Zr/Fe] 


[Ba/Fe] 


[Nd/Fe] 


23439A 
23439R 
+5°3640 


-1.05 
-1.11 
-1.27 


+0.30 
+0.39 
+0.53 


+0.28 
+0.35 
+0.38 


+0.38 
+0.57 
+0.59 


+0.23 
+0.20 
+0.38 


+0.24 
+0.39 



1 



X 
3 



> 
'■+-' 

a 

(U 

ct: 




Cas 
Teff 5350 
[Fe/H] -0.92 



HD 65583 
Teff 5305 
[Fe/H] -0.82 



HD 108564 
Teff 4660 
[Fe/H] -1.18 



7799 



[Rb/Fe] 
0.0, 0.3, 0.6 




[Rb/Fe] 
0.0, 0.3, 0.6 




[Rb/Fe] 
0.0, 0.3, 0.6 



Si I Rb 



7800 
Wavelength (A) 



7801 



Fig. 8.— ACTUAL FIGURE 1 
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Fig. 10.— ACTUAL FIGURE 3 




Fig. 11.— ACTUAL FIGURE 4 



T 



X 
3 



> 
'■+-' 

a 

(U 

ct: 




1.1 




.8 



HD 103095 
Teff 5040 
[Fe/H] -1.48 



HD 23439A 
Teff 5200 
[Fe/H] -1.05 



HD 23439B 
Teff 4950 
[Fe/H] -1.11 




6134 



Zr I 



V I 



6135 
Wavelength (A) 



6136 



Fig. 13.— ACTUAL FIGURE 6 
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[Fe/H] 



Fig. 14.— ACTUAL FIGURE 7 



